Spin order and lattice frustration in optimally doped manganites. A high temperature 

NMR study. 
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The physics underlying the complex glassy phenomena, which accompany the formation of po- 
larons in optimally doped manganites (ODM) is a cumbersome issue with many unexplained as- 
pects. In this article we present ^''^La and ^^Mn NMR in the temperature range 80K - 900K of 
ODM Lao.erCao.aaMnOs. We show that local lattice distortions, established in the Paramagnetic 
(PM) phase for T < 700K, induce a genuine spin-glass state, which for T < Tc consolidates with 
the Ferromagnetic (FM) state into a single thermodynamic phase. Comparative NMR experiments 
on Lao.77Cao.23Mn03, Lao.59Cao.4iMn03, and Lao.roSro.aoMnOa demonstrate the dominant role of 
lattice distortions, which appear to control (i) the stability of the spin glass phase component and 
(ii) the kind (1st or 2nd order) of the PM-FM phase transition. The experimental results are in 
agreement with the predictions of the compressible random bond - random field Ising model, where 
consideration of a strain field induced by lattice distortions, is shown to invoke at a discontinuous 
(1st order like) change of both the FM and the "glassy" Edwards- Anderson (EA) order parameters. 

PACS numbers: 75.47.Lx, 75.50.Lk, 76.60.-k, 75.30.Et 
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I. INTRODUCTION 

The study of strong electron correlations in transi- 
tion metal oxides unveiled a complex world of interweav- 
ing properties, concerning their spin, charge, and crys- 
tal structure. Predominant examples arc high tempera- 
ture superconducting cupratcs and hole doped mangan- 
ites. Competition among different interactions in these 
systems generates spectacular phenomena, such as the 

mesoscopic 



formation of charge and spin stripe: 
phase separation^ii, and the colossal magnetoresistance 
(OMR) effeclj^i^. At the same time frustration of in- 
teractions gives rise to the appearance of freezing and 
glassinessiSiiij expressed with slow relaxation, aging, and 
other signatures of glassy systemsi^ii^iii^. However, it is 
not yet clear, whether this kind of glassiness is a funda- 
mental property of strongly correlated electron systems, 
or the consequence of quenched disorder, which produces 
uncorrelated charge and spin fluctuations. 

In the case of hole doped manganites, exemplified by 
the prototype Lai_a;Ca2;Mn03 (LCMO) family, substitu- 
tion of La'^^ ions with the divalent alkaline-earth metal 
ion Ca^+. invokes the replacement of Jahn- Teller (JT) 
active Mn'^+ ions with JT inactive Mn'*+ ions. Here, 
frustration is generated by competition between (i) co- 
herent Jahn- Teller (JT) lattice distortions, which favour 
charge localization, and (ii) the Double Exchange mech- 
anism, which favours motion of Cg electrons between 
adjacent FM ordered Mn'^^ and Mn'*+ ions, and tends 
to smooth out lattice distortions. In the doping range 
(0.2 < a; < 0.5) the FM metallic phase is imposed, and 
the ground state is that of an isotropic FM metal (Fig- 
ure [1]). However, a number of experiment a^^i^^'^^ have 
shown that at optimally doping, x « 0.33, in the region 
of Tc, nanoscalc lattice polarons (local lattice and spin 
textures) are formed, which drive the PM to FM phase 
transition to 1st order— i^i^^. At the same time, neu- 



tron scattering experiments, have shown that polarons 
might form a glassy phase for T > Tc, which subse- 
quently melts to a polaron fluid by further increasing 
temperaturei^i^ (Figure [T|). Other experiments revealed 
that the anomalous 1st order "like" phase transition is ac- 
companied by frequency and temperature dependences^, 
and strong relaxation effects2ii2^, which are reminisnece 
of relaxor ferroelectrics and spin glasses. Similar effects 
have been observed in Ndo.yoSro.soMnO^^. By contrast, 
no such effect have been observed in other ODM man- 
ganites, such as Lao.ToSro.soMnOs (LSMO(0.30)), which 
exhibit a conventional 2nd order phase transition^Si^iSl. 

It is evident by now that the remarkable differences be- 
tween the Ca-doped and Sr-doped systems are imposed 
by their different high temperature crystal symmetry: In 
LCMO(0.33), a Rhombohedral to Orthorombic (R-3c to 
Pnma) phase transition occurs by coohng, at « 
accompanied by the onset of strong collective JT dis- 
tortions, while LSMO(0.33) remains in the Rhombohe- 
dral R-3c phase at all temperatures^^ (inset of Figure 
[T]), where static JT displacements are forbidden by the 
crystal symmetry. However, it is unclear whether the 
observed "glassiness" in LCMO(0.33) is simply due to 
random freezing of polaronic distortions, or a collective 
spin-glass transition of the Ising type, driven by frus- 
trated interactions and internal stresses. 

Another puzzling issue is that by shifting away from 
X = 0.33, where the stronger glassy phenomena are ob- 
served, the 1st order PM-FM phase transition turns to 
2nd order (Figure[T|). Figure [2] demonstrates H/M vs 
isotherms of LCMO(0.23, 0.33, 0.41) and LSMO(0.30) in 
the vicinity of Tc. It is clearly seen that the panel for 
X = 0.33 shows negative slope in the lower Af^ region, 
which according to the Banerjee criterio n^°i^° is a clear 
sign that LCMO(0.33) belongs to the 1st order transition. 
By moving away from optimal doping, in both directions, 
the slope becomes positive signifying that the transition 
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LSMO(0.30) systems. 



FIG. 1: The phase diagram of Lai-jjCaj^MnOg for 0.0 < 
X < 0.5. The polaron glass and dynamic polaron regimes 
are defined according to ref.— . The green points (line) de- 
fine the PM to FM transition line in magnetic field 9.4 
Tesla. The inset shows the corresponding phase diagram for 
Lai-ajSr^^MnOs from ref.—. The blue line is the transition 
line between the R-3c and Pnma crystal structures. 



becomes of the 2nd order. However, in the lower dop- 
ing regime, scaling arguments, which characterize 2nd 
order phase transitions are not directly applicable, and 
there are questions about the exact nature of the phase 
transitional. 

In order to shed light to this difficult issue, we have per- 
formed ^^^La and ^^Mn NMR in the temperature range 
80 K to 900 K, on optimally doped LCMO(G.33), as well 
as LCMO(0.23), LCMO(0.41), and LSMO(0.30) powder 
samples. In contrast to diffraction techniques, where long 
range order with a coherence length of at least lOOnm is 
required, NMR is a local probe and therefore ideal for 
the characterization of short range ordered systems. As 
far as we know, this is the first time that NMR experi- 
ments on manganites were performed at temperatures as 
high as 900K. Our experiments suggest that in case of 
LCMO(0.33) strong local lattice distortions, apparently 
of polaronic origini^, are formed below the Rhombohe- 
dral to Orthorombic transition temperature at w TOOK. 
These distortions are responsible for the appearance of a 
spin-glass state in the PM phase, while for T < Tc the 
spin glass state merges with the FM phase, comprising 
a new thermodynamic phase (a kind of collective spin- 
glass) . The order in this novel phase resembles the order 
of spin-glasses in a magnetic field, i.e. the long range FM 
order and the short range glassy order coexist. Contrary 
to reentrant spin-glass phases, by lowering temperature 



the growing FM order is shown to truncate the glassy 
phase component. 

This intriguing phase transition is excellently de- 
scribed by the compressible random interaction-random 
field Ising model^, where strains induced by polaronic 
distortions are taken into consideration. In case of 
LCMO(0.33) a discontinuous change of the spin and EA 
order parameters at are foreseen by the model, while 
for T < Tc the EA order parameter qea decays by de- 
creasing temperature, in agreement with the experimen- 
tal results. We notice that the EA order parameter qEA 
in the FM state is defined hy qEA = <1ea — M'^, where 
QEA = jfJ2i{Sif = [{Si) ]av The spin operator Si 
takes the values ±1, (...) represents the time average, 
and [....]a-u denotes the disorder average. Experiments 
show that by moving away from optimal doping, lattice 
distortions faint out and the glassy component weakens 
rapidly, while the absence of an internal strain field is re- 
sponsible for turning the phase transition from 1st order 
to 2nd order, as predicted by the model. 



II. MATERIALS AND METHODS 

Lai-ajDiajMnOs (Di=Ca, Sr) samples were prepared by 
thoroughly mixing high purity stoichiometric amounts of 
CaCOa (SrCOg), LaaOg, and MnOa. The mixed pow- 
ders formed in pastile form, reacted in air at 1400 °C 
for several days with intermediate grinding and then 
slowly cooled down to room temperature. X-ray diffrac- 
tion measurements were performed on a D500 SIEMENS 
diffractometer, showing that all samples are single phase 
materials with very good crystallinity. Magnetization 
measurements were performed on a Quantum Design 
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MPMSR2 superconducting quantum interference device 
magnetometer. At temperatures higher than 350K, mag- 
netization measurements were performed on a Lakeshore 
3700 vibrating-sample magnetometer (VSM). The NMR 
experiments were performed on a home built broadband 
spectrometer operating in the frequency range 5 — 800 
MHz, in 9.4 Tcsla and zero external magnetic fields. 
Spectra were acquired by the spin-echo point by point 
method while varying the frequency, because of the 
large spectral width of the resonance lines. An Oxford 
1200CF continuous flow cryostat was employed for mea- 
surements in the temperature range 80 — 350 K and an 
Oxford HTIOOOV furnace for measurements in the range 
300 - 900 K. 

III. HIGH TEMPERATURE La NMR 

The local magnetic and structural properties of the 
PM phase for all samples have been investigated by ap- 
plying ^'^^La NMR in magnetic field 9.4 Tesla, at tem- 
peratures as high as 900K. In the presence of an ex- 
ternal magnetic field B, ^^^La(/ = 7/2) nuclei expe- 
rience the Zeeman interaction, which splits the degen- 
erate nuclear energy levels into 2/ + 1 equidistant en- 
ergy levels, with energies Em = m^HB. In addition, 
the ^^^La nucleus is coupled to the local electric field 
gradient (EFG) tensor, through its electric quadrupole 
moment Q. We notice that a site with cubic symme- 
try has an EFG equal to zero; so the size of the effect 
of the electric quadrupole interaction on the NMR spec- 
trum measures the degree of deviation from cubic sym- 
metry of the surroundings of the nuclear site. In the 
presence of the Zeeman and Quadrupolar interactions the 
frequency of the NMR spectrum corresponding to transi- 
tions between levels m and to — 1 is given by, LUm ~ "fB + 
l inf{2?-i) (3cos26i - 1 - risin^ecos'^4>) (2m - 1). Here, e 
is the charge of the electron, eq is equal to V^z, r] is the 
asymmetry parameter of the EFG tensor 77 = ^ 
and 6 , 4> are the angles between the principal axis of the 
EFG tensor and the magnetic field. In case of LCMO 
manganites, the magnetic field B at the site of the La 
nuclei is equal to the sum of the external field B and the 
transferred hyperfine field Bhf = {l/jfijA < S >, where 
A is the hyperfine coupling constant and < S > the aver- 
age electronic spin of the eight nearest Mn neighbors^^iS. 
In the PM phase, large frequency shifts proportional to 
the magnetic susceptibility are produced by the hyper- 
fine coupling in the presence of an external magnetic field. 
According to the above equation, the transition between 
the TO = +1/2 and to = —1/2 levels (the so called cen- 
tral transition) is unaffected by the electric quadrupole 
interaction to 1st order, while the distance in frequency 
of all other transitions (sattelites) from the central tran- 
sition, depends solely on the quadrupolar interaction. In 
a powder sample, the crystal axes, and hence the EFG, 
are distributed at random angles with respect to the ap- 
plied magnetic field. The angular average of the satellite 
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FIG. 3: i^^La NMR spectra for optimally doped LCMO(0. 33), 
in 9.4 Tesla external magnetic field, as a function of temper- 
ature. Spectra are presented relatively to the central line 
frequency uq. The dotted line shows schematically the bor- 
der line between the R-3c and Pnma phase regimes. The left 
inset depicts the calculated powder pattern spectrum in com- 
parison to the experimental spectrum at T = 730K. The right 
inset shows a comparison of spectra at 80K (open circles) and 
400K (solid line). 



patterns gives rise to a characteristic frequency distribu- 
tion. The size of this distribution is proportional to the 
EFG (and consequently to the lattice distortions) and 
independent of the magnetic field. 

Figure [3] demonstrates ^'^^La NMR spectra of the 
LCMO(0.33) system in the temperature range 320 — 
900K. At temperatures higher than 700K spectra exhibit 
typical powder patterns, consisting of the narrow cen- 
tral line (sa 150kHz broad), and the broad satellite fre- 
quency distribution (SFD) (« 1.6MHz broad); the lat- 
ter is the fingerprint of the EFG in the i? — 3c Rhom- 
bohedral phase. The inset at the top left side of Fig- 
ure [3] depicts the calculated powder pattern spectrum 
at 730K, which matches excellently with the experimen- 
tal spectrum. For the simulation an electric quadrupo- 
lar coupling vq — = 21MHz, asymmetry param- 
eter n = 0.5 and a dipole-dipole interaction of the La 
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nuclear sites with the Mn ions of w 20kHz were con- 
sidered. By decreasing temperature, at T « TOOK, a 
significant increase of the SFD width is observed, which 
marks the transition from the i? — 3c to the Pnma crys- 
tal structure^, and the appearence of strong incoher- 
ent JT displacements. This structural phase transfor- 
mation is shown by synchrotron radiation measurements 
to be independent of the applied mag nctic field2^. At- 
tempts to simulate the powder pattern of the sattelites 
for T < TOOK were only possible by assuming a distri- 
bution of EFG tensor components, i.e. a distribution 
of lattice distortions. Most remarkably, by further de- 
creasing temperature, the average width of the SFD re- 
mains invariant down to the lowest measured tempera- 
ture, as clearly seen in the inset at the top right side of 
Figure [31 Hence, from the viewpoint of the La site, the 
size and distribution of lattice distortions in the Pnma 
phase does not change by decreasing temperature. At the 
same time, the narrow central future starts to broaden 
by cooling and disappears at « 320K, in agreement with 
previous measurements^. The extreme broadening of 
the central line, which up to 1st order depends solely 
to magnetic interactions, is a clear sign about the onset 
of a broad distribution of magnetic susceptibilities in the 
Pnma phase (the NMR signal broadening is proportional 
to the product of the magnetization and the suscepti- 
bility distribution width). Evidently, the appearance of 
short range correlated polaronsi^, produces a distribu- 
tion of exchange couplings, which is directly reflected on 
the spin order of the system. 

Figure S] shows ^^^La NMR spectra of the 
LCMO(0.23), LCMO(0.41), and LSMO(0.30) sys- 
tems in the temperature range 300 — 885K. In case 
of LCMO(0.23), at high temperatures spectra depict 
the typical NMR powder pattern as in LCMO(0.33). 
However, two basic differences between the two systems 
are that (i) the SFD of LCMO(0.23) is at aU measured 
temperatures sufficiently narrower than the SFD of 
LCMO(0.33), and (n) the SFD width of LCMO(0.23) 
varies smoothly by decreasing temperature in the PM 
phase, while in LCMO(0.33), after increasing abruptly at 
T w TOOK, it remains invariant by lowering temperature 
(inset in Figure [5]). At first sight this is a surprising 
result, because LCMO(0.23) exhibits (i) transition from 
the Pnma to the R — 3c phase at approximately the 
same temperature as LCMO(0.33), and (ii) neutron 
scattering PDF measurements have shown that local JT 
distortions decrease by increasing Ca doping^. A pos- 
sible explanation is that close to doping x = 0.25, small 
polarons might form an ordered lattice, where Mn'*+ 
ions have 6 neighbouring JT distorted Mn"'"'^ sites, fitting 
together in a space filling 3D network, which minimizes 
lattice distortions and strains^. The picture of reduced 
lattice distortions in the Pnma phase of LCMO(0.23) is 
supported by the caculated powder pattern spectra. The 
right inset in Figure S] depicts the calculated spectrum 
in comparison to the experimental one at 410K. For 
the fit a quadrupolar coupling constant lyq — 28MHz, 
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FIG. 4: i^^La NMR spectra for LCMO, x = 0.23(left), 0.41 
(middle), and LSMO x = 0.3 (right) at various temperatures. 
Spectra are presented relatively to the central line frequency 
vo- The left inset shows a comparison of LCMO x — 0.23 
(blue), and 0.33 (black) experimental spectra at T = 410K, 
and 400K respectively. The right inset depicts the calculated 
powder pattern spectrum in comparison to the experimental 
spectrum for x = 0.23 at T = 410K. 



an asymmetry parameter rj = 0.6 and a dipolc-dipole 
interaction of the La nuclear sites with the Mn ions 
w 50kHz were considered. In addition, a large Gaussian 
broadening was necessary for the fit, which is indicative 
of the presence of substantial lattice distortions. Nev- 
ertheless, these distortions are sufficiently weaker than 
in LCMO(0.33), where in the Pnma phase adequate fit 
is only possible by considering overlapping spectra with 
a broad distribution of quadrupolar coupling constants. 
On the other hand, in case of LCMO(0.41) the SFD is 
almost absent and only the magnetic broadening of the 
central line is observed by decreasing temperature. This 
result is in agreement with the phase diagram in Figure 
[TJ which shows no static polarons for x = 0.41, but 
rather the presence of dynamic polaronsi^ii^ with a zero 
time average of local lattice distortions in the time scale 
of NMR experiments. Similar spectra but for different 
reasons are shown by LSMO(0.30). This system remains 
at all temperatures in the i? — 3c Rhombohcdral phase 
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FIG. 5: The ^'^''La NMR frequency as a function of tempera- 
ture for LCMO(0.23) (•), LCMO(0.33) (•), LCMO(0.41) (o), 
and LSMO(0.30) (•). In the inset the width of the NMR 
frequency distribution is shown for (i) LCMO(0.23) central 
line (•), and satellite powder pattern (o), (ii) LCMO(0.33) 
central line (•), and satellite powder pattern (o), and (iii) for 
LSMO(0.30) central line (•), and satellite powder pattern (o). 




240 280 320 360 400 440 

Frequency [MHz] 



I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' 

240 280 320 360 400 440 

Frequency [MHz] 



^^Mn NMR spectra (open circles) of LCMO(0.33) 
and 200K and LSMO(0.30) (left figures) at 80K 



FIG. 6: 

at 80K 

and 300K, respectively. The spin-spin T2 relaxation times 
(crosses) as a function of frequency are also shown in the same 
Figures. By increasing temperature spectra shift to lower 
frequencies, while for LCMO(0.33) they become asymmetric 
broader on approaching Tc from below. The T2 minimum at 
the peak of the spectra at 80K is due to the Shul-Nakamura 
inetractions. 



(inset in Figure [T|) , where static JT distortions are by 
symmetry forbidden, and therefore it depicts only a very 
narrow SFD, while the central line broadens gradually 
by lowering temperature. 

Figure shows the ^^^La NMR frequency at 9.4 Tcsla 
as a function of temperature, for all four measured sys- 
tems. The critical temperatures Tc, defined as the in- 
flection point on the ly vs. T curves, are found to 
be 260K for LCMO(0.23), 300K for LCMO(0.33) and 
LCMO(0.41) and 400K for LSMO(0.30). These values 
(shown as green points in Figure [T]) define the PM to 
FM transition line at 9.4 Tesla in the T — x phase di- 
agram, which is parallel shifted to higher temperatures 
from the transition line in zero external magnetic field. 
We notice that although in 2nd order phase transitions 
critical temperature does not exist in the presence of a 
magnetic field, the observed temperature variation of the 
NMR frequency is fairly steep, which is reminicence of a 
true phase transition. Furthermore, it is observed that 
the transfered hyperfine field at saturation, expressed by 
the ^•^^La NMR frequency at 80K, depends almost lin- 
early on the doping concentration, in agreement with 
previous ^'^^La measurements in zero external magnetic 
field^. The inset in Figure [5] shows the full width at half 
maximum (FWHM) of the central line (filled circles) and 
the FWHM of the SFD (open circles) for LCMO(0.23), 
LCMO(0.33), and LSMO(0.30). As previously discussed, 
in case of LCMO(0.33) the transition from the i? — 3c to 
the Pnma phase by cooling, is clearly seen to occur at 



T « TOOK. Below that temperature, the abrupt increase 
of the SFD width declares the appearance of strong static 
local JT distortions, which remain invariant down to the 
lowest measured temperature. It is also clearly seen that 
the SFD width and thus the local lattice distortions, are 
in LCMO(0.33) by far the largest among the four mea- 
sured systems. Another important notification is that 
contrary to LCMO(0.33) where the central line disap- 
pears at T « 300K, apparently due to extreme inhomo- 
geneous broadening, in the rest three systems the central 
line remains visible at all temperatures, and its width 
varies with temperature in exactly the same way as the 
^^^La NMR frequency. This is a clear sign that - unlike to 
the LCMO(0.33) system - there is no broad distribution 
of magnetic susceptibilities. 

IV. ^^Mn NMR MEASUREMENTS IN ZERO 
EXTERNAL MAGNETIC FIELD 

The invisible with ^^^La NMR "magnetic" glassy phase 
component for LCMO(0.33), is possible to monitor by 
applying ^^Mn NMR in zero external magnetic field. For 
reasons of comparison NMR data are compared with the 
prototype (in respect to the 2"*^ order nature of the phase 
transition) LSMO(0.30) system in the temperature range 
80K up to temperatures close to Tc- 

^^Mn NMR in zero external magnetic field probes di- 
rectly the electron spin state of single Mn ions through 
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Mocssbauer spectroscopy from ref<^ are also presented in 
the same figure. An excellent match between the NMR 
and the Moessbauer hyperfine fields is observed. An- 
other important experimental feature is that in case of 
LCMO(0.33) Bhf decreases abruptly to zero at Tc, a clear 
mark of the 1st order nature of the phase transition. On 
the contrary, in case of LSMO(0.30) the hyperfine field 
decreases continously to zero by approaching Tc from be- 
low, as expected for a 2nd order phase transition. 

V. THE COMPRESSIBLE 
RANDOM-INTERACTIONS AND RANDOM 
FIELDS ISING MODEL 



FIG. 7: The hyperfine filed Bhf as a function of temperature, 
obtained from '"^^Mn NMR measurements, for LCMO(0.33) 
(o) and LSMO(0.30) (o). The filled circles (•, •) are the 
corresponding Bhf values as obtained by Moessbauer spec- 
troscopy on lightly Sn-doped samples from ref.— . The solid 
lines are theoretical fits by applying the self consistent equa- 
tions (2) and (3). The inset shows the corresponding spectra 
for LCMO(0.33) (blue lines) and LSMO(0.30) (black lines). 



the hyperfine field B^f ~ {l/jh)A < S >, and therefore 
it is possible to resolve the different Mn charge states, 
i.e., localized Mn'^+, Mn'^+ and the intermediate FM va- 
lence states. Figure ini shows ^^Mn NMR spectra for both 
samples at 80K and 200K (BOOK), respectively. Spectra, 
as expected, were found to consist of a single line that is 
caused by motional narrowing, due to the fast electron 
hoping between the Mn'^"'', Mn^+ manganese ions with 
frequency much higher than the NMR frequency^!. In 
the same figure the nuclear spin-spin relaxation time T2 
as function of frequency is demonstrated. T2 in man- 
ganites depends strongly on temperature^^, apparently 
due to fluctuations of the hyperfine field caused by hop- 
ing of the electron holes. At lower temperatures, an 
additional contribution to T2 is observed, due to the 
Suhl-Nakamura relaxation mechanism^S, which produces 
a characteristic T2 minimum at the center of the spec- 
tra, in agreement with previous results^. This fre- 
quency dependence of T2 disappears by increasing tem- 
perature. In case of LSMO(0.30) T2 is shown to be 
frequency independent on approaching Tc from below. 
However, in LCMO(0.33) another frequency dependent 
spin-dynamics appears, which is expressed as a mono- 
tonic increase of T2 by increasing frequency. This effect 
indicates the presence of a distribution of low frequency 
fluctuations close to Tc in LCMO(0.33), which is absent 
in LSMO(0.30). At the same time, for T > 150K a sig- 
nificant asymmetric broadening of the lincshape is ob- 
served in LCMO(0.33), which is absent in LSMO(0.30), 
in agreement with previous NMR works^^iiS. 

Figure [7| demonstrates lineshapes and the normalized 
hyperfine field Bhf/Bhf{0) for both LCMO(0.33) and 
LSMO(0.30) as a function of temperature. For rea- 
sons of comparison the Bhf / Bhf{Q) values obtained with 



The discontinuous hyperfine field variation at Tc and 
the inhomogeneous ^^Mn NMR lineshape broadening for 
T < Tc observed in LCMO(0.33), are nicely explained 
by considering an Ising Hamiltonian, where random- 
interactions and random-fields, together with a strain 
field induced by strong local lattice distortions are taken 
into consideration^li^. 



J^j + ^ Jlj e u 



SiSj — ^ ^ fiSi 

i 

+ (1) 



The last term expresses the macroscopic elastic energy. 
The random interactions Jf^ and random fields fi are 
assumed to be independently distributed according to 
Gaussian probability densities. This is an extension of 
the Shcrrington-Kirkpatrick model^, which in the weak- 
disorder limit predicts coexistence of spin-glass and FM 
order.^^. 

In the replica-symmetric mean field approximation^, 
the reduced magnetization and the EA order parame- 
ters are possible to be calculated by the self-consistent 
equation 
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where, Tg, Ta and T^ are temperatures, which char- 
acterize the variance of the random interactions and 
random fields distributions, and the strength of the 
strain field, respectivcly^i. By solving these equations 
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FIG. 8: Experimental ^^Mn NMR lineshapes (left) and the- 
oretical simulations (right) for LCMO x = 0.33 at various 
temperatures, as obtained by formula (4). 



self-consistently, it is possible to fit the experimen- 
tal Bhf I Bhfo vs. T curves in Figure [71 In case of 
LCMO(0.33) data are fitted by considering = 265K, 
Tg = 40K, Ta = 5K, and Te = 75K. This solution has 
a jump at Tc, in agreement with the experimental data. 
On the contrary, LSMO(0.30) is nicely fitted with a 2nd 
order phase transition, by ignoring both randomness and 

strain field. For T > Tc eq. 3 implies that qea = i-C- 
the EA order parameter appears to decay by increasing 
temperature. This is in agreement with the ^^^La NMR 
results presented in Figure [3l where the appearance and 
gradual narrowing of the central-transition NMR line at 
temperatures higher than 320K, indicates the supression 
of the spin-glass phase component at elevated tempera- 
tures. In a similar way, the effective EA order parameter 
qea is shown to decay rapidly in the FM phase by de- 
creasing temperature (Figure [7]) . 



On the basis of the fitting parameters, it is posible to 
calculate the temperature variation of the frequency dis- 
tribution, /(^), which is proportional to the probability 
distribution of the local magnetization VF(m}^, 



/(i/) cx W{'m) 
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rp2 



2^2 



1/2 I 



■ exp 



[arctanh{m) - - 



Qea- 



T2 



(4) 



By using eq. 4, and the fitting parameters from Figure 
[7] an excellent agreement is obtained between the exper- 
imental and the simulated LCMO(0.33) lineshapes, as 
shown in Figure [51 Evidently, the ^^Mn NMR hneshape 
broadening observed in LCMO(0.33) by approaching 
from below, is a sign of coexistence of long range FM 
order and short range spin-glass order, which accompa- 
nies the appearance of polaronic lattice distortions in this 
temperature region. 



VI. CONCLUSIONS 

On the basis of our experimental results we anticipate 
that optimally doped LCMO(0.33) shows a unique mag- 
netic behaviour in comparison to all other measured sys- 



tems, by varying the temperature across Tc. Specifically, 
by entering the Orthorombic Pnma crystal structure at 
T K, 700K on cooling, a broad distribution of strong lat- 
tice distortions appears, apparently of polaronic origini^, 
which in turn induces an equally broad inhomogeneous 
distribution of Mn electron spin polarizations. Short 
range correlations of such polaronic distortionsi^ii^ ap- 
pear to give rise, to a spin-glass state, which for T < Tc 
coexists with the FM order, but it gradually disappears 
by further decreasing temperature. We emphasize that 
this genuine "collective" spin-glass state is not produced 
by quenched disorder, but is self generated and depends 
only on the competition among generic interactions such 
as magnetic exchange and JT interactions. Important 
role is played by the magnitude of lattice frustration, 
which creates an inhomogeneous strain field that controls 
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the stability of the glassy state and the kind of the mag- 
netic phase transition (1st or 2nd order). It is interesting 
to note that in the simplified model Ising Hamiltonian we 
have used, simulations show that the induced randomness 
on effective J-couplings dominates in comparison to ran- 



dom fields created by trapped polaronic charge carriers. 
Indeed, changes of the J couplings (produced by local 
variations of the Mn — O bond lengths) are expected to 
have a sufficiently stronger effect than the weak random 
dipolar fields. 
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